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Abstract 

We have developed a series expansion method for calculating the uniform 
magnetization, M, as a function of the applied field h for Heisenberg systems 
at T = 0. The method involves introducing Ising anisotropy along the z-axis 
together with an applied uniform field along the x-axis. On extrapolation to 
the isotropic limit, one recovers the magnetization for the Heisenberg system 
with an applied field along the x-axis. This method circumvents the difficulties 
in developing perturbation theory associated with the commuting nature of 
the uniform field. The results developed for two-chain ladders appropriate for 
the material Cu2(C5Hi2N2)2Ci4 are in good agreement with the experimental 
data. In addition, uniform susceptibility is calculated by high temperature 
expansions and also compared with the experimental data. 

PACS Indices: 75.10.-b., 75.10J., 75.40.Gb 
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Recently several spin-ladder materials have been discovered, which exhibit properties 
characteristic of one-dimensional systems The material (^(CsH^^^C^ represents 

a particularly simple Heisenberg spin-system from a theoretical point of view as it forms || a 
spin-ladder with the inter-chain coupling, J± much larger than the intra-chain coupling, J , 
making it ideally suited for perturbation theory and for testing out quantitative theoretical 
calculations. For this system, the nuclear relaxation rate 1/T\ along with the Knight shift 
and field dependent magnetization have been measured by Levy et al. Here, we would like 
to study this system by series expansion method and in particular try out a new method for 
calculating field dependent properties of the Heisenberg spin systems. 

In recent years series expansion methods have played a significant role in studying mag- 
netic properties of quantum spin models and real materials f4j||. In addition to high temper- 
ature expansions, a variety of T = perturbation expansions have been developed to study 
quantum spin Hamiltonians. Important technical advances have been made in calculating 
the excitation spectra by these methods ||. 

One general problem that, to our knowledge, has not been considered before is related to 
dependence on properties associated with a commuting variable. Example of such variables 
are chemical potential and uniform magnetic field. The difficulty posed by such variables 
is that at T = different values of magnetization or density lead to different sectors of 
the Hilbert space, which are not connected to each other. Thus, naively, one would need 
to consider all different sectors of the Hilbert space to find the global minimum and hence 
determine the dependence on these fields. This is clearly a formidable task. 

Here we circumvent this difficulty by explicitly breaking the symmetry at intermediate 
stages of the calculation. For the Heisenberg-Ising model in a transverse field, the field 
term no longer commutes with the Hamiltonian. By series extrapolation, we recover the 
properties of the Heisenberg limit, where the field term does commute with the Hamiltonian. 
It is assumed that by the limiting procedure, implicit in the series analysis, the global 
ground-state with respect to all possible magnetizations is obtained in the Heisenberg limit. 
Comparison with experimental data confirms that this is a valid procedure for calculating 
the field dependent properties. 

We begin with the spin-ladder Hamiltonian in an external field 

H = H l + h-'£S i (1) 

i 

where Hi represents the spin-ladder Hamiltonian in zero field. High temperature expansions 
for the system with Jj_ = J.,, as well as T = Ising and Dimer expansions for this model 
for h — have been discussed earlier |7|. 

Here, in order to compare with the experimental data for (^(CsH^^^CU, and deter- 
mine the exchange constants for the system, we calculate the high temperature series for 
uniform susceptibility x(T) with arbitrary y 2 = J../J±. The series coefficients are given in 
table 1. To make the comparison easier, we take J±/J,, = 5, a value inferred earlier by 
Levy et al. ||. Since the susceptibility is presented in arbitrary units, we are left with two 
fitting parameters: an overall magnitude of x is se t by a free parameter, whereas the fit 
of the temperature scale is controlled by the exchange constant J±. We find that the best 
fit is obtained with Jj_ = 13.96K. Comparison with experimental data is shown in Fig. 1 
and verifies previous assertions that the above model with J_l/J,, = 5 is appropriate for the 
material under consideration. 
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In this paper we are primarily interested in T = field dependent magnetization. The 
properties in zero field including the spin-gap is most accurately calculated by the dimer 
expansions. The gap is estimated to be A = 4. 10855 J.,. However, as discussed earlier, this 
expansion is not suitable for calculating the field dependent properties. Instead, we consider 
the following model which depends on the parameter A, 

H = H Q + XV (2) 

1=2 

Ho = J a ^M^M+i + ^-L ^liS^i + tY(-l) l+l Si i 

i,l=l i l,i 

1=2 

T/ T \^ I qx qx , qy qy \ , j V~~V qx qx , qy qy \ ± sr^( 1 Q z \ h Q x 

V — J a \ J l,i°l,i+l "r °l,i J l,i+l) "T J -L Z^V°l,i°2,j D l,i°2,i) l H l ) J l,i T a / , 

i,l=l i 2,i l,i 

where S^j denotes the S=l/2 spin at the ith site of the Zth chain. Notice that the staggered 
field terms (the terms with coefficient t) are simply added to improve convergence, and 
we adjust the coefficient t to get the smoothest terms in series with a typical value being 
t = 2Jj_ . For A = 1 the model reduces to the spin-ladder Heisenberg model in a uniform 
field. However, for A < 1, the field term does not commute with the Hamiltonian and thus 
mixes states with different magnetization. In practice, the only complication this brings is 
that the dimensionality of Hilbert space for a given graph that needs to be considered is 
increased. 

We have obtained the expansions for uniform magnetization M, the expectation value 
of l/NJ2i,i Sfi i n the ground state, up to order A 14 for a range of field values (the series 
coefficients for several h/J,, values are presented in Table 2), the calculations involved a 
list of 2767 linked cluster of up to 14 sites. The resulting series are extrapolated to the 
Heisenberg point by Pade and differential approximants ||. The resulting plots for the 
magnetization are shown in Fig. 2. It is encouraging that the magnetization appears to go 
to zero at the onset field, h c \ = 4.108557^ known accurately from the gap calculated from 
the dimer expansions. Although, the dependence of the magnetization on the field is linear 
for most of the field values, near the onset field h&, it can be fit to a power law behavior 

M = c[{h - h cl )/j/. (3) 

Using the accurately known h c \ and performing two-point fits leads to the estimates f3 = 
0.55(10) and c = 0.18(3), consistent with theoretical expectations 0. Near the upper critical 
field h C 2, where the magnetization saturates to the maximum value of 1/2, our results become 
less accurate, making it difficult to obtain any critical exponents. The ratio h C 2/h c \ is ~ 1.8 
in our calculations, which is close to the experimental value of 1.7. The full plots are in 
reasonable agreement with the experimental data. 

Conclusions: In this paper, we have developed a series expansion method for calculating 
the field dependent properties of the Heisenberg model. The method is shown to work well for 
the spin-ladder system and provides a good description for the material Cu2(C5Hi 2 N 2 ) 2 Cl4. 
In future, we hope to use this method to study the field dependent properties of the Heisen- 
berg models in other geometries. 
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FIGURES 



FIG. 1. Comparison of the calculated susceptibility \ f° r the system with J // /J± = 0.2 (shown 
as the solid lines) with the experimental data (shown as open symbols) of Levy for the material 
Cu 2 (C 5 H 12 N 2 ) 2 Cl 4 . 

FIG. 2. Uniform magnetization M for the system with J,,/J± = 0.2 as a function of 
applied field h/J /f obtained by Ising expansion, the solid line is the power law behavior 
M = 0.18[(/i — h c i)/J ,] - 55 estimated by the two-point fits to the region near h c \. 
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TABLES 



TABLE I. Series coefficients for high-temperature expansion of the uniform susceptibility 
Tx{T) = J2n,m a n,m(y2) n l3 m of the system with 7/2 = J // /J±, and (3 = J±/T. Nonzero coeffi- 
cients a ntm up to order m = 9 are listed. 



(n,m) 


dn,m 




(n,m) 






(n,m) 


dn.m 




(n,m) 






( 0, 0) 


2.500000000x10" 




( 0, 5) 


2.115885417x10" 




( 1, 7) 


2.556694878x10" 




( 7, 


8 ) 


-3.694080171x10" 




(0, 1) 


-6.250000000x10" 


2 


( 1, 5) 


-1.790364583x10" 


-3 


(2, 7) 


-5.249023438x10" 


-4 


( 8, 


8) 


2.766200474x10" 


-4 


(1,1) 


-1.250000000x10" 


1 


( 2, 5) 


3.255208333x10" 


-3 


(3, 7) 


6.212022569x10" 


-4 


( 0, 


9 ) 


3.899804709x10" 


-6 


(0, 2) 


-1.562500000x10" 


2 


( 3, 5) 


-2.604166667x10" 


-3 


(4, 7) 


-4.855685764x10" 


-4 


( 1, 


■>) 


-2.440800742x10" 


-5 


(1.2) 


6.250000000x10" 


2 


(4, 5) 


9.765625000x10" 


-4 


( 5, 7) 


5.045572917x10" 


-4 


( 2, 


9) 


4.863436260x10" 


-5 


(0, 3) 


1.302083333x10" 


3 


( 5, 5) 


-1.367187500x10" 


-3 


( 6, 7) 


9.087456597x10" 


-5 


( 3, 


9) 


-6.667041274x10" 


-5 


(1,3) 


7.812500000x10" 


3 


( 0, 6) 


-1.044379340x10" 


-4 


( 7, 7) 


6.200396825x10" 


-5 


( 4, 


9 ) 


1.052977547x10" 


-4 


(2, 3) 


-1.562500000x10" 


2 


( 1, 6) 


2.766927083x10" 


-4 


( 0, 8) 


1.513768756x10" 


-6 


( 5, 


9) 


-1.893664163x10" 


-4 


(3, 3) 


1.041666667x10" 


2 


( 2, 6) 


-3.255208333x10" 


-5 


( !, 8) 


1.850430928x10" 


-5 


( 6, 


9 ) 


1.067227157x10" 


-4 


(0, 4) 


1.627604167x10" 


3 


( 3, 6) 


1.736111111x10" 


-4 


( 2, 8) 


-1.121399895x10" 


-4 


( 7, 9) 


-6.055075025x10" 




( 1,4) 


-5.208333333x10" 


3 


( 4, 6) 


-1.155598958x10" 


-3 


( 3, 8) 


1.617916047x10" 


-4 


( 8, 


9 ) 


-1.049223400x10" 


-4 


(2, 4) 


3.906250000x10" 


3 


( 5, 6) 


1.497395833x10" 


-3 


( 4, 8) 


5.907331194x10" 


-5 


( 9, 




5.028403415x10" 


-5 


(3, 4) 


-5.208333333x10" 


3 


( 6, 6) 


-1.082356771x10" 


-3 


( 5, 8) 


-2.031598772x10" 


-4 










(4, 4) 


3.255208333x10" 


3 


(0, 7) 


-3.986661396x10" 


-5 


( 6, 8) 


4.343184214x10" 


-4 











TABLE II. Series coefficients of the uniform magnetization M = J2 i a^A 1 of the system with 
J///J± = 0-2 and staggered field t = 2J±. Nonzero coefficients a\ up to order i = 14 are listed. 



i 


a,i for h 1 J a = 4 




di for h/J a = 4.5 


for hj J a = 5 




di for h/J a = 5.5 


1 


-2.35294117647x10" 


1 


-2.64705882353x10" 


1 


-2.94117647059x10- 


1 


-3.23529411765x10" 


1 


2 


-1.09573241061x10" 


2 


-1.23269896194x10" 


-2 


-1.36966551326x10- 


2 


-1.50663206459x10" 


2 


3 


9.14025678774x10" 


2 


1.06396227962x10" 


-1 


1.22649301848x10- 


1 


1.40301724400x10" 


1 


4 


6.20129611723x10" 


2 


5.97803768408x10- 


-2 


5.40239556347x10- 


2 


4.43521601235x10" 


2 


5 


1.13232334931x10" 


2 


-6.87114677881x10- 


-3 


-3.16482548625x10- 


2 


-6.35515201898x10" 


2 


6 


1.51540112742x10" 


3 


-9.22803967212x10- 


•3 


-2.09238359642x10- 


2 


-3.20631922064x10" 


2 


7 


1.04579690591x10" 


2 


1.06631743925x10" 


-2 


1.41218986085x10- 


2 


2.30445918289x10" 


2 


8 


7.27178133729x10" 


3 


4.84071702453x10- 


-3 


1.81299605913x10- 


3 


-2.22929845852x10" 


3 


9 


-2.65316666928x10" 


3 


-9.18394611547x10- 


-3 


-1.77224406625x10- 


2 


-2.87388172848x10" 


2 


10 


-3.02918907022x10" 


3 


-3.20313568381x10- 


-3 


4.56086296539x10" 


4 


1.06088934260x10" 


2 


11 


4.16113147230x10" 


3 


1.07015220405x10" 


-2 


2.17700230213x10- 


2 


3.73144477607x10" 


2 


12 


7.17721455928x10" 


3 


9.49330584821x10- 


-3 


7.35902021044x10- 


3 


-5.77787480240x10" 


3 


13 


3.76830635331x10" 


3 


-1.96870341084x10" 


-3 


-1.59885104470x10- 


2 


-4.04664282941x10" 


2 


14 


4.80548328343x10- 


4 


-5.66412721220x10- 


3 


-1.13503013822x10- 


2 


-6.54745511421x10" 


3 


i 


m for h 1 J a = 6 




di for h/J a = 6.5 


di for h/ J a = 7 




di for h/ J a = 7.5 


1 


-3.52941176471x10- 


1 


-3.82352941176x10" 


-1 


-4.11764705882x10- 


1 


-4.41176470588x10" 


1 


2 


-1.64359861592x10" 


2 


-1.78056516724x10" 


-2 


-1.91753171857x10- 


2 


-2.05449826990x10" 


2 


3 


1.59493430485x10" 


1 


1.80364354971x10" 


-1 


2.03054432723x10- 


1 


2.27703598609x10" 


1 


4 


3.03734528767x10- 


2 


1.16962964638x10" 


-2 


-1.20708465458x10- 


2 


-4.13195135825x10" 


2 


5 


-1.03065008581x10" 


1 


-1.50607485884x10" 


-1 


-2.06526481612x10- 


1 


-2.71092352584x10" 


1 


6 


-4.06274351241x10- 


2 


-4.40370259369x10- 


-2 


-3.91005725854x10- 


2 


-2.19638419584x10" 


2 


7 


4.00968938376x10- 


2 


6.83808087160x10- 


-2 


1.11397495384x10- 


1 


1.72990002765x10" 


1 


8 


-8.31253081902x10- 


3 


-1.83810619203x10- 


-2 


-3.56778569062x10- 


2 


-6.52026110830x10" 


2 


9 


-4.29821115807x10- 


2 


-6.15590391942x10- 


-2 


-8.59943598221x10- 


2 


-1.18256118000x10" 


1 


10 


3.08422130226x10- 


2 


6.59751313197x10- 


-2 


1.22493770522x10- 


1 


2.09158221300x10" 


1 


11 


5.57220816744x10- 


2 


7.29635570762x10- 


-2 


8.15164933911x10- 


2 


6.90272394747x10" 


2 


12 


-4.00161047870x10- 


2 


-1.09477296280x10- 


-1 


-2.32383122267x10- 


1 


-4.30523552978x10" 


1 


13 


-7.38343538295x10" 


2 


-1.06305578450x10- 


-1 


-1.13164071654x10- 


1 


-4.53616541698x10" 


2 


14 


2.82846389387x10- 


2 


1.25033040054x10- 


-1 


3.29080311091x10- 


1 


6.96568120714x10" 


1 
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